Finding the elusive sliding phase in superfluid-normal phase transition smeared by c-axis disorder 
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We consider a system composed of a stack of weakly Josephson coupled superfluid layers with c-axis disorder 
in the form of random superfluid stiffnesses and vortex fugacities in each layer as well as random inter-layer 
coupling strengths. In the absence of disorder this system has a 3D XY type superfluid-normal phase transition 
as a function of temperature. We develop a functional renormalization group to treat the effects of disorder, and 
demonstrate that the disorder results in the smearing of the superfluid normal phase transition via the formation 
of a Griffiths phase. Remarkably, in the Griffiths phase, the emergent power-law distribution of the inter-layer 
couplings gives rise to sliding Griffiths superfluid, with an anisotropic critical current, and with a finite stiffness 
in a-b direction along the layers, and a vanishing stiffness perpendicular to it. 



The interplay of disorder and broken symmetry remains a 
challenging and relevant problem for correlated quantum sys- 
tems. The effects of disorder in one dimensions, where the 
effects of quantum fluctuations are enhanced, is most dra- 
matic, giving rise to Anderson localization [1], Dyson singu- 
larities, and random singlet phases [2]. Recent studies, both 
experimental and theoretical, concentrated on the superfluid- 
insulator transition of Bosonic chains [3], and strongly argued 
that disorder alters the universality of that transition [4, 5]. 
While uncorrected disorder in higher dimensions has a less- 
ened effect, we must raise the question: how does correlated 
disorder, which only varies in a subset of directions, affects 
thermal and quantum phase transitions in higher dimensions? 

In this work, we study this question by concentrating on the 
superfluid insulator transition in 3D Bose gases, that is split 
into a series of pancake clouds by a ID optical lattice with 
disorder which varies only along the lattice direction, but not 
parallel to the clouds. While this question is of much theo- 
retical interest, and is now also of experimental relevance as 
we outline below, it was not addressed so far. The effects of 
the disorder could be as mundane as just shifting the transi- 
tion point, or as important as resulting in a new universality 
class of the transition or obliterating it altogether. Indeed, we 
shall show that the interplay between disorder along the c- 
axis and the a-b plane Berezinskii-Kosterlitz-Thouless (BKT) 
physics [6, 7], smears the transition giving rise to an interme- 
diate Griffiths phase [8, 9] that occupies a wide region of the 
phase diagram. Furthermore, in a subphase within this Gri- 
iffith phase, the superfluid becomes split into an array of 2D 
puddles that have no phase coherence along the c-axis, thus 
realizing the illusive sliding phase paradigm [10], supporting 
superflow only in the a- and b- but not c-directions. 

The questions we raise are fast becoming important for ex- 
periments. Experiments on ultracold atoms observed both the 
BKT transition in large 2D "pancakes" produced by very deep 
ID optical lattices [11], and Anderson localization of Bosons 
in ID disordered optical lattices [3, 12]. The system we study 
here can be realized by constructing a stack of large 2D "pan- 
cakes" using a disordered ID optical lattice and tests the ef- 
fects of disorder near the 2D-3D crossover [13, 14]. 

The model which we analyze and describe below consists 




FIG. 1: Left: Schematic diagram of the model: red ovals (purple 
bars) represent the superfluid layers (Josephson couplings) with size 
inversely proportional to vortex fugacity £ m (directly proportional to 
Josephson coupling J m ). The effects of the real space RG is to merge 
strong layers, and decimate weak layers. Emergent length-scale Lj 
corresponds to the typical separation between strong layers. 
Right: Schematic diagram of the RG flows showing the superfluid 
(SF) and normal fixed points along with the Griffiths fixed line. The 
black dashed line represents physical configurations, with points to 
the right corresponding to higher temperatures. The Griffiths fixed 
line is split into two segments, corresponding to the regimes with fi- 
nite and zero c-axis superfluid response. The star indicates a possible 
unstable fixed point [20]. 



of a set of coupled 2D superfluid layers. Each layer has a su- 
perfluid stiffness K m , vortex fugacity (akin to vortex density 
per coherence length) ( m , and Josephson coupling (to the next 
layer) J m . K m , £ m and J m are initially random and uncorre- 
cted [21], see Fig. la. To analyze this model, we combine a 
Kosterlitz-Thouless like momentum space renormalization for 
the in-plane degrees of freedom [7, 15] with a real-space RG 
(see, e.g., [2, 5]). In the real space-RG decimation, strongly 
coupled layers ( J m ^ 1) are merged, while vortex-ridden lay- 
ers (£ m ~ 1) are considered to be essentially normal and are 
perturbatively eliminated. 

Before plunging into the analysis, let us summarize the 
phase diagram we find, see Fig. lb and Table I. At low temper- 
atures the system forms a 3D superfluid. As the temperature 
is raised, a Griffiths phase appears; in it, the system breaks up 
into 2D superfluid puddles, each composed of one or several 
"pancakes", with weak (power law distributed) inter-puddle 
tunneling. As the temperature is increased further, the c-axis 
superfluid response disappears altogether, while the system 
remains superfluid in the a- and b-directions, realizing a slid- 
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TABLE I: Phase diagram indicating the properties of the various 
phases. 



ing phase. At yet higher temperatures, the in-plane superfluid 
response smoothly vanishes as the system becomes fully nor- 
mal. 

The Griffiths phase is perhaps the most surprising aspect of 
our results. At intermediate temperatures, the flow leads to 
a fixed line characterized by a stationary power-law distribu- 
tion of inter-layer tunnelings, J m , which are P(J) ~ J" 7-1 
(Fig. lb). The appearance of these Griffiths phase power laws 
are a direct consequence of the disorder. Most layers have 
strong fluctuations, and turn insulating; Neighboring layer of 
either side, can still exchange bosons but with a smaller am- 
plitude, e.g., Jeff = Jm-i ■ Jm if layer to is eliminated. The 
elimination of all incoherent layers marks a first epoch in the 
RG flow, and upon its end the internally coherent layers are 
separated from each other by, on average, Lj incoherent lay- 
ers Fig. la. Lj determines vy. vj ~ log[l/J]/ij, with J 
the typical initial Josephson coupling. In the subsequent RG 
epoch, layers only merge, but vj remains unchanged. 

The Griffiths phase can be separated into two regimes. 
A sliding regime with vj flowing to vj < 1 (as indicated 
in Fig. lb), where there is no c-axis stiffness, and a Grif- 
fiths superfluid with a finite c-axis stiffness and vj > 1. 
Both regimes have a vanishingly small c-axis critical cur- 
rent. To wit, the critical current of n layers is determined 
by the weakest effective tunneling between them. The ex- 
pectation value for the longest run of weak layers is lZ n ~ 
log 1/jw [n(l - p„eak)] [16] (with p weak the probability of a 
layer to be normal in the first epoch; Lj = (1 — p„ ea k) 1 3> 

/ \ Lj log J 

1). The weakest link is thus I c ~ ( J 

Model — Let us now describe the model and its anal- 
ysis more precisely, before discussing more of its conse- 
quences and experimental implications. Following Ref. [15], 
our model consists of a set of coupled 1 + 1 dimensional 
(Euclidean) sine-Gordon models with partition function Z = 
Tr exp - J2 m (S s G;m + <Sj ; m,m+i) where 
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S S G;m is the sine-Gordon action that describes the density 
waves and vortices in the m-th layer; <Sj ;m . m +i is the to 
to to + 1 tunneling; 9 m (x,y) and 4> m (x : y) are the super- 
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FIG. 2: Schematic representation of relevance of J's and C's as a 
function of temperature. 



fluid order-parameter phase variable and its conjugate, respec- 
tively ([0 m (r), d x <p m >(r')] = iir5(r - r')S m . m >). We define 
J m = Jm/T, K m = K. m /T, and Cm ~ exp(-£' core , m /T) 
as the reduced Josephson coupling, superfluid stiffness, and 
vortex fugacity at temperature T, where E cole m is the vortex 
core energy. We note that to define the sine-Gordon model 
we must specify the short-distance cut-off scale. We choose a 
single cut-off a for all layers, and work in the units in which 
a = 1. 

Renormalization Group — Our analysis relies on a com- 
bined c-axis real space and a-b momentum space RG. The 
momentum space transformation is given by [22]: 
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(4) 
(5) 



To lowest order in Cm and J m , there is a range of super- 
fluid stiffnesses 1/4tt < K m < 2/7r in which both the vor- 
tex fugacity Cm and the Josephson coupling J m are relevant. 
The competition between the two gives rise to the Griffiths 
phase. Outside this range the system is either strongly super- 
fluid (large K m ) or strongly insulating (small K m ), Fig. 2. 

As the in-plane momentum shell RG proceeds, the real- 
space RG (RSRG) merges layers where J m rises to 1, or elim- 
inates layers where vortex fugacity Cm reaches to 1. When 
a Josephson coupling becomes large, J, = 1, the relative 
phase AO = Qi+\ — B L of the two neighboring layers becomes 
locked and the two layers merge into a single super-layer hav- 
ing K eS = Ki + K i+1 and Ceff = Q ' d+i- Similarly, if 
one of the vortex fugacities becomes large, Q — 1, then the 
conjugate field <fii in that layer becomes locked and vortices 
proliferate. Upon integrating out the incoherent layer, we find 
that it suppresses tunneling across it to J e ff = Jj_i • Jj. These 
RG rules make it convenient to parametrize J and C in terms 
of their logs j = log(l/J), z — log(l/C) which yields: 



Jeff = +ji 



Z e ff 



-1 + Zi 



(6) 



Next, instead of a numerical analysis of the RG outlined 
above (which we will fully pursue in a separate work [17]), 
let us use the RG procedure to derive the approximate flow of 
the distribution functions for K, j, and z, and their universal 
aspects. First, note that the RSRG layer merging step leads to 
strong correlations of K's and C's. Therefore, alongside the 
distribution Pj for jm's, we use the joint probability distri- 
bution Q Z K for z m 's and AT m 's. The fRG equations resulting 
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from Eqs. (5, 6) and the: 
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4^ 4xK 2 y i-yj stands 
for g®(g) with being the step function; dl and Qi are short- 
hand for dK\dz\ and Q z ^ i where it is unambiguous. Briefly, 
the first terms of Eqs. (7) and (8) correspond to the action of 
the linear in J m and ( m terms of the momentum space RG 
Eqs. (3) and (4). The remaining terms correspond to the ac- 
tion of the real space RG, where we keep in mind the fact that 
the distributions must be normalized. The normalization is 
accomplished by rescaling the distributions Pj and Q Z K when- 
ever layers are removed from the system. The fRG equations 
must be supplemented by absorbing wall boundary conditions 
Q K = and P Q - = which remove the small z's and j's 
(large £'s and J's) from the distributions when layers are dec- 
imated or merged. In order to compute physical observables 
we also keep track of n(£), the number of surviving layers at 
RG scale I: 



dn 
~dt 



= ~n ^hP + J dK x {2 - irK^Q ^ 



(9) 



Note that the structure of the fRG and of the resulting flows 
are similar to those in Ref. [9, 18] for the damped transverse 
field Ising model [23]. 

To study the evolution of Pj(l) and Q Z K {1) under coarse 
graining, we numerically integrate Eqs. (7) and (8). To 
parametrize the initial distributions at temperature T (and 
length scale a), we choose smooth functions with the follow- 
ing bounds: 0.04 < TJ < 0.11, 1.0 < TK < 1.5, and 



.67T-1.5/T < Q < 



-1.67T-1.0/T 



[24]. 



Results — A numerical analysis of the flows reveals three 
phases: (1) superfluid phase - all layers merge, (2) Griffiths 
phase - power law distributions Pj <~ J VJ ~ X with a finite v j, 
(3) insulating phase - all layers decimate. Phases (1) and (3) 
correspond to the usual superfluid and insulating fixed points, 
the Griffiths phase, regime (2), corresponds to a new fixed line 
that is induced by disorder [20]. 

Within the Griffiths phase, the flow of Pj and Q Z K distri- 
butions occurs in two epochs, as depicted in Fig. 3. In the 
first epoch both mergers and layer eliminations take place, 
which quickly results in the formation of power law distri- 
butions with flow of all three exponents. However, as merg- 
ers lead to ever increasing K's while eliminations do not, the 



flow eventually exhausts all weak layers by eliminating them, 
and only the strongly superfluid layers remain. In the second 
epoch, only J's (but not £'s) remain relevant as all the sur- 
viving K's exceed 2/ir, so while layers continue to merge, no 
eliminations occur. As a result, i/j saturates while both vr 
and i/£ decay to zero exponentially in the fRG flow parame- 
ter vk ~ ~ e~ 2vjl - We see, therefore, that the Griffiths 
fixed line corresponds to a line of fixed i/j's with vk — > and 
-> 0. 

Both the in-plane p a b and the out-of-plane p c superfluid 
responses have a very peculiar behavior within the Griffith 
phase, and can be used to as probes. The mean values of the 
superfluid responses may be obtained from the distributions 



via 



Pab = n j dzdKKQ z K , (10) 
p c = x (l) + e- 21 (nj dje j P^j , (11) 



where, n, which is given by Eq. (9), is the fraction of surviving 
layers and is needed to normalize the response to the scale of 
the original system. x(7) obeys 



dt X (!) = nP e 2l h 



x(o) = 0, 



(12) 



and accounts for stiffness within the superfluid "puddles". As 
p c depends on the area, we include the factor of e~ 2 ' in Eq. 1 1 
to account for its renormalization. We plot p a b and p c as a 
function of T near saturation (at large value of the RG pa- 
rameter £) in Fig. 4. The smearing of the phase transition is 
reflected in p ab decreasing smoothly, as the temperature is in- 
creased, until reaching zero at the end of the Griffiths fixed 
line (without following any power law). On the other hand, 
p c decreases much faster, becoming zero within the Griffith 
phase at the point where lim^oo vj{t) becomes smaller than 
unity. The disappearance of p c signals the onset of the elusive 
sliding subphase of the Griffith phase, see Table I. 

Experimental probes — In our analysis we found that dis- 
order smears the superfluid-to-normal phase transition. This 
could be probed experimentally in the ultra-cold atom setting. 
The superfluid response as well as the critical current could 
be measured by jolting the confined gas (e.g. quickly displac- 
ing the trap potential) and looking at the decay of the cen- 
ter of mass oscillations [19]. Alternatively, one could look at 
correlations by removing the optical lattice and the trap po- 
tentials and allowing the atoms from the various pancakes to 
expand and interfere. The key signature of the Griffiths phase, 
in interference experiments, is very strong shot noise which 
results from the interference of several weakly coupled super- 
fluid droplets [17]. Alternatively, in the mesoscopic setting, 
the Griffiths phase could appear in artificially grown struc- 
tures composed of alternating layers of superconducting and 
insulating films of varying thicknesses. In this setting the 
anisotropy superfluid responses and critical currents could be 
measured directly. 
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FIG. 3: Left three panels: flow of exponents vq, vk, and vj under the action of the coarse graining transformation for three different initial 
distributions corresponding to temperatures T = 3.0, 3.2, and 3.4. I ~ 3 separates the first epoch, in which all three exponents flow from the 
second epoch in which only and vk flow. The asymptotic value of vj at long length-scales indicates that T = 3.0 corresponds to a Griffiths 
superfluid, while T = 3.2 and 3.4 correspond to Griffith insulator. 

Right panel: semi-log plot of typical distributions (from top to bottom) Q[z], R[K], and P\j], obtained by solving Eqs. (7) and (8) numerically 
and the corresponding exponential fits (black solid lines) that are used to obtain the values of exponents i>k, and vj. The inset depicts the 
distribution R[K] on a linear-linear plot, the oscillations arise from additions of K values when layers merge. 
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FIG. 4: The in-plane (p a b) and out-of-plane (p c ) stiffness as a func- 
tions of temperature evaluated using value of RG parameter is suf- 
ficiently large (I = 10) so that the flows are saturated. While p a b 
undergoes a smeared phase transitions, as indicated by the absence 
of a critical point with power law behavior, p c undergoes a continu- 
ous phase transition at T ~ 3.1. 

Conclusions — In this manuscript we investigated the ef- 
fects of reduced dimensionality disorder on a phase transition 
at a higher dimensionality. We focus on a model that we be- 
lieve is relevant to experiments in ultra-cold gases in optical 
lattices, and mesoscopic systems such as stacked supercon- 
ducting films. Naively, one would expect that, the strong dis- 
order picture of Ref. [2] should be in effect. However, the 
classification scheme of Ref. [8] indicates that the existence 
of the BKT transition for a single 2D layer boosts the impor- 
tance of disorder in our system, resulting in the stronger effect 
of the smearing of the phase transition. Using a functional 
renormalization group scheme that we develop, we show that 
this is indeed the case. Further, we show that in the transi- 
tion region the system becomes essentially two dimensional. 
We find that the reduction of dimensionality is reflected in the 
strong anisotropy of physical observables like critical current 
and superfluid response. 

As we were finalizing this manuscript, we became aware of 
a complementary investigation of the random superfluid stack 
from a scaling perspective by Vojta and Narayanan, which is 
consistent with our findings. 
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